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I.  Introduction 

The  DMSP  F7  spacecraft  was  launched  on  18  Nor,  1983  into  a  98.74  degree 
inclination  orbit,  with  apogee  844  kn  altitude  and  perigee  822  kn.  (Rich, 
1984) .  The  primary  purpose  of  the  spacecraft  was  to  obtain  tropospheric 
meteorological  data.  However,  a  triaxial  flux gate  magnetometer  was  included 
on  the  spacecraft  in  order  to  monitor  the  geophysical  environment.  This 
report  deals  with  the  examination  of  these  magnetometer  data  to  evaluate 
their  usefulness  in  describing  the  earth’s  core-produced  geomagnetic  field. 

The  field  is  assumed  to  be  curl  free  and  so  representable  by  a 
potential  function  in  the  form  of  a  spherical  harmonic  series  (Langel,  1987) 

n*  n 

1)  V  a  a  E  E  (a/rju+lfgn^osm^  +  hnmsinm^]Pn®(cos5) 

n=l  m=0 
n**  n 

+  a  E  E  (r/a)&[qn®cosm^  +  s^sinm^JPnOCcosfl) 
n=l  m=0 

where  a  is  the  mean  radius  of  the  earth  (taken  to  be  6371.2km),  r,  0,  and  l 
are  the  standard  spherical  coordinates,  and  the  Pn»(cos0)  are  the  Schmidt 
quasi-normalized  form  of  associated  Legendre  functions  of  degree  n  and  order 
a.  The  summation  involving  (a/r)  represents  sources  internal  to  the 
satellite  orbit,  whereas  that  involving  (r/a)  represents  external  sources 
such  as  the  ring  current. 

The  magnetic  field  is  then  given  by 

2)  B  =  -TV 

Theoretically,  1)  and  2)  hold  only  if  n*  and  n**  go  to  infinity.  The 
measured  internal  B  contains  contributions  from  both  the  Barth’s  core  and 
from  its  crust;  n*  in  1)  is  thus  chosen  so  that  V  represents  fields  from  the 
core  but  not  the  crust.  Langel  and  Estes  (1982)  concluded  that  the  core 
field  dominates  for  n<13  and  the  crustal  field  for  n>15.  The  n*  chosen  for 
highly  accurate  measurements,  such  as  from  the  MAGSAT  satellite,  usually 
equals  13,  but  because  the  DMSP  observations  are  less  accurate,  n*  does  not 
exceed  11  in  any  DMSP  model. 

The  DMSP  F7  magnetometer  was  mounted  on  the  satellite  body,  as  opposed 
to  being  attached  to  a  boom,  because  of  spacecraft  engineering  constraints 
(Rich,  1984).  The  magnetometer,  a  triaxial  fluxgate,  was  aligned  wi  ;h  the 
spacecraft  X,  Y  and  Z  axes,  which  are  defined  as  follows:  X  is  vertically 
down,  Y  is  along- track  and  Z  is  cross- track.  The  three  sensor  units  were 
built  by  the  Schonstedt  Instrument  Co.,  Reston,  VA,  in  the  1960s.  The 
electronics  unit  for  the  magnetometer  was  built  by  the  Applied  Physics 
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Laboratory  of  Johns  Hopkins  University,  Laurel,  HD,  based  on  the  design  of 
the  MAGSAT  fluxgate  aagnetoaeter . 

The  aagnetoaeter  acquired  field  aeasureaents  at  a  rate  of  20  saaples 
per  second.  Measureaents  were  in  the  fora  of  counts,  with  one  count 
equalling  12  nano-Teslas  (nT) .  According  to  Rich  (1984) ,  the  instruaent  was 
not  intended  to  survey  the  aain  geomagnetic  field,  so  it  was  not  calibrated 
with  high  accuracy  on  the  ground,  nor  recalibrated  in  orbit. 

Because  of  the  close  proxiaity  of  the  aagnetoaeter  to  on-board 
electronic  instruaentation,  its  data  were  contaminated  by  non-random 
instrumental  noise,  with  magnitudes  of  up  to  several  thousand  nT.  The 
attitude  of  the  spacecraft  was  measured  to  an  accuracy  of  about  0.1  degree, 
or  360  arc-seconds.  While  this  attitude  accuracy  is  not  as  good  as  that 
obtained  with  MAGSAT,  in  principle  it  is  of  sufficient  accuracy  to  enable 
meaningful  vector  measurements.  In  the  absence  of  other  near-Earth 
satellite  magnetic  field  data  for  this  time  period,  and  in  view  of  the 
success  of  methods  used  on  MAGSAT  to  solve  for  spacecraft  fields  (see 
section  III) ,  it  was  decided  to  investigate  the  possibility  of  processing 
the  DMSF  F7  data  to  a  stage  where  they  aay  be  useful  for  main  field 
modeling. 

The  remainder  of  this  report  will  deal  with  several  topics  related  to 
the  DMSP  data:  1)  the  processing  steps  and  problems  encountered  in  removing 
both  randoa  and  coherent  noise  from  the  data,  2)  the  correction  of  data 
offsets  and  attitude  errors  using  the  GSFG  main  field  software  (FIT 
program),  3)  the  utilization  of  processed  data  in  modeling  the  core-produced 
geomagnetic  field,  and  4)  the  evaluation  of  those  models  in  regards  to  the 
usefulness  of  this  type  of  data  for  main  field  studies. 


II.  Transformation  of  On- tape  Data  to  Magnetic  Readings. 

The  magnetoaeter  data  contained  on  the  basic  DMSP  data  tapes  received 
from  the  Air  Force  is  in  the  fora  of  magnetometer  counts,  which  must  be 
converted  into  field  values  in  nT  in  order  to  be  useful.  Data  is  arranged 
on  the  tape  as  a  header,  containing  time  and  position  information  for  each 
minute  of  operation,  followed  by  60  magnetometer  readings  (1  per  second) . 

All  times  on  the  data  are  rounded  to  the  nearest  second.  All  positions  on 
the  original  Air  Force  data  are  expressed  in  nautical  miles,  rounded  down  to 
the  nearest  nautical  mile.  Nautical  miles  were  converted  to  kilometers  (one 
nautical  mile  equals  about  1.8  km)  in  subsequent  data  processing. 

Rounding  of  the  time  and  position  introduces  an  error  of  unknown 
magnitude  into  the  data.  An  estimate  of  the  possible  magnitude  of  this 
error  can  be  made  as  follows.  The  satellite  velocity  is  about  7.4 
km/second.  A  one  second  time  error  is  therefore  equivalent  to  an  along- 
track  position  error  of  7.4  km.  Combined  with  the  possible  1.8  km  roundoff 
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error  in  position,  the  total  along- track  position  error  may  be  as  such  as 
0.2  km.  These  errors  are  in  addition  to  those  due  to  actual  error  in  the 
computed  epheaerides. 


Lange 1  (1978)  has  estiaated  the  effects  of  orbit  error  on  near-earth, 
polar,  satellite  aagnetic  field  experiaents.  His  results  are  suamarized  in 
Table  1: 


Table  1 

Component  Maximum  gradient  (nT/ka) 

Vertical  Along- trk  Cross-trk 


Mean  of  gradient  magnitude  (nT/ka) 
Vertical  Along-trk  Cross-trk 


Br 

-28.0 

-13.3 

8.8 

14.9 

5.3 

2.0 

B<? 

18.0 

-8.5 

23.4 

7.8 

3.0 

1.7 

8.4 

-2.0 

23.3 

2.6 

0.7 

2.1 

B 

-28.0 

-8.1 

5.7 

18.4 

2.5 

1.6 

A  position  error  of  1.8  ka  in  the  vertical  direction  will  give  maximum 
field  errors  of  50.4,  32.4,  and  15.12  nT  in  Br,  3# ,  and  Bi,  respectively.  A 
position  error  of  1.8  ka  in  the  cross-track  direction  will  give  maximum 
errors  of  12.24,  42.12,  and  42  nT,  and  an  error  of  9.2  km  along-track  will 
give  maxiaua  field  errors  of  123,  60,  and  19  nT  in  Br,  B$,  and  B^, 
respectively.  It  seeas  probable  that  the  resulting  errors  are  in  the  50-150 
nT  range  and  that  they  may  be  at  least  partly  systematic  in  nature. 

20  readings  per  second  for  each  of  the  magnetometer  X,  Y  and  Z  axes 
were  originally  recorded  by  DMSP,  although  only  readings  #1  and  #11  were 
written  on  the  tapes  sent  to  the  Geology  and  Geomagnetism  branch  at  Goddard 
(See  Appendix  A  for  on- tape  data  formats) .  For  the  Goddard  main  field 
studies,  only  the  magnetometer  reading  associated  with  the  header  record  was 
utilized.  This  preserved  sufficient  data  density  (one  reading  per  minute) 
to  fully  describe  the  main  geomagnetic  field. 
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The  magnetometer  was  calibrated  prior  to  launch  at  the  NASA  Goddard 
Space  Flight  Center  aagnetic  test  chamber  with  the  following  results  (Rich, 
1984): 


3)  Measurement  = 

(nT) 

Calibration  Matrix  * 
(nT/count) 

Measurement 

(counts) 

+  Bias 
(nT) 

(radial)  fBxl 

1  i 

ri2.1001  -0.0055  0.01931 

|  I 

rcxi 

1  1 

r-0. 06531 

I  1 

1  1 

(along- trk)  IByl  = 

1  1 

1-0.0247  12.1883  -0.0101 1 
|  | 

1  1 

*  ICy  1 

1  1 

+  159.87331 
|  | 

1  1 

(cross-trk)  IB*J 

(.  0.0069  0.0232  12.1735J 

ICjjJ 

139.4228J 

This  equation  is  used  to  compute  the  magnetic  field  in  spacecraft 
coordinates  in  nT,  given  a  reading  in  magnetometer  counts.  However,  this 
calibration  does  not  take  into  account  the  field  from  the  spacecraft,  which 
adds  greatly  to  the  bias  vector.  This  vector  must  be  determined  from  in¬ 
flight  data.  Later  work  using  the  FIT  program  (see  Field  Value  corrections) 
accomplished  this,  and  re-determined  the  bias  vector  as:  (89nT,  8457nT, 
-1441nT)  for  radial,  along- track  and  cross-track  measurements,  respectively. 
These  values  are  still  somewhat  approximate  and  require  small  corrections 
which  will  be  introduced  and  discussed  in  sections  VIII  and  IX.  In 
addition,  for  computational  ease,  it  was  decided  to  redefine  the  spacecraft 
system  to  be  compatible  with  the  MAGSAT  coordinate  system,  so  that  the 
spacecraft  X  axis  is  defined  as  cross-track,  Y  is  radially  down  and  Z  is 
along- track.  The  correct  transformation  of  DMSP  magnetometer  counts  to 
nano-teslas  in  spacecraft  coordinates  compatible  with  MAGSAT  is  thus: 


4)  Measurement  =  Calibration  Matrix  *  Measurement  +■  Bias 


(nT) 

(nT/ count) 

(counts) 

(nT) 

(cross-trk)  fBxl 

I  | 

fO.0069  0.0232  12.1735  1 

I  1 

rcxi 

1  1 

r-i44ii 

1  1 

(radial)  IByl 

1  1 

=  112.1001  -0.0055  0.0193  1 

|  | 

1  1 

♦  1  Cy  1  + 

1  1 

1  89  1 

1  | 

1  1 

(along-trk)  IB*) 

1-0.0247  12.1883  -0.0101J 

IC*I 

l  8457J 

All  data  discussed  in  the  remainder  of  this  report  are 
processed  through  this  equation  and  have  units  of  nT. 
measurement  vector  in  equation  #4  will  henceforth  have 
meaning  the  vector  is  in  MAGSAT  spacecraft  coordinates 
through  the  corrections  discussed  in  section  III. 


assumed  to  have  been 
The  (Bx,By,Bz) 
the  label  Bspu, 
and  not  yet  processed 
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III.  Field  Value  Corrections. 

According  to  Rich  (1984),  the  DUSP  aagnetoaeter  may  be  misaligned 
relative  to  the  spacecraft  by  as  much  as  0.5  degree  per  axis,  with  the 
misalignment  measured  to  an  accuracy  of  about  0.1  degree.  Also,  bending  of 
the  spacecraft  body  may  result  in  further  misalignment.  In  addition,  the 
values  of  the  three  magnetic  components  may  be  in  error  by  a  fixed  bias  or 
by  a  multiplying  factor.  The  FIT  program  has  the  capability  to  solve  for 
corrections  in  these  parameters  in  conjunction  with  the  least  squares  main 
field  solution.  The  theory  of  this  adjustment  is  as  follows  (see  also 
Estes ,  1983) : 

The  FIT  program  computes  three  types  of  adjustments  to  vector 
satellite  magnetometer  data:  1)  A  diagonal  calibration  matrix  containing 
aslopea  parameters,  which  is  multiplied  times  the  measured  vector  to  correct 
for  magnetometer  drift,  2)  a  bias  correction  vector  which  is  subtracted  from 
the  measured  vector  to  correct  for  constant  magnitude  offsets,  and  3)  a 
rotation  matrix  which  is  multiplied  times  the  measured  vector  to  correct  for 
angular  offsets  of  the  aagnetoaeter  from  ideal  satellite  coordinates.  These 
adjustment  parameters  are  applied  to  the  measured  uncorrected  data  in 
spacecraft  coordinates  according  to  the  equation: 

5)  Bspc  =  TSR  *  TCAL  *  (Bspu  -  bias) 

where:  Bspu  is  the  uncorrected  measurement  vector  in  spacecraft 
coordinates ,  as  given  in  equation  4) . 

Bspc  is  the  corrected  measurement  vector  in  spacecraft 
coordinates . 

bias  is  a  vector  of  magnetometer  bias  corrections  in  addition 
to  those  given  in  equation  4) . 

TCAL  is  the  calibration  correction  matrix  of  slope  parameters. 

TSM  is  the  rotation  correction  matrix. 


The  elements  of  bias  are :  (BSi ,  BS2 ,  BS3) ,  where  BSi  are  component 
biases  derived  in  the  FIT  program,  with  values  derived  and  discussed  in 
sections  VIII  and  IX. 


TCAL  has  elements: 

fl/SLi 

1 

0  0  1 
1 

!  0 

1 

1/SL2  0  1 
1 

l  0 

0  I/SL3J 

where  SLi  are  slopes  derived  in  the  FIT  program.  SLi  and  BSi  are 
applied  to  the  satellite  X  axis,  SL2  and  BS2  to  satellite  Y  axis,  and  SL3 
and  BS3  to  the  satellite  Z  axis  components. 
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The  elements  of  TSMl  ere  based  on  three  Euler  angles  (  ex»Cy>  and  eB) 
solved  in  execution  of  the  FIT  program.  (Note:  In  the  FIT  program,  as  of 
2/28/88,  ex  is  denoted  e2,  ey  is  denoted  ei,  and  eB  is  denoted  £3.)  Using 
the  notation  TSMj.  j ,  where  i  is  the  matrix  row  and  j  the  matrix  column,  these 
are: 

6)  TSMn  =  cosey *cose  B 

TSMl 2  =  cos£y*cos£x*sin€B  +  sinCy*sinex 
TSMl 3  =  -cosey*sinex*sineB  +  siney*cosex 
TSM21  =  -sinew 

TSM22  =  cosex+cosez 

TSM23  =  -sin£X*C08£w 

TSM3I  a  -sin£y*COS£jj 

TSM32  =  cos6y*sin6x  _  siney*sineB*cosex 

TSM33  a  cos£y*cos£x  +  siney*sinex*sineB 

The  TSM  matrix  in  the  FIT  program  is  derived  from  3  rotation  matrices 
(denoted  Rx,  Ry,  RB)  in  the  spacecraft  coordinate  system.  Rx  is  a  left- 
handed  rotation  through  the  angle  £x,  about  the  spacecraft  X  axis.  Ry  is  a 
left-handed  rotation  about  the  spacecraft  Y  axis,  with  angle  of  rotation  £y. 
Rs  is  a  right-handed  rotation  through  the  angle  eB,  about  the  spacecraft  Z 
axis. 


The  matrices  Rx>  Ry  and  RB  are  thus 


7) 

rx  =  ri 

1 

10 

1 

0 

0  1 

1 

-sinex 1 

1 

Ry  =  fcosCy 

1 

0  siney  1 

1 

COS£x 

1  0 

I 

1  0  1 

1 

10 

sinex 

cosexJ 

1 

L-siney 

1 

0  COSCyJ 

Rb  =  f  cos£b  sin£B  0  1 

I  I 

I -sine B  coseB  0  I 

I  I 

10  0  1  J 


Rotations  provided  by  Rx,  Ry  and  RB  are  illustrated  in  Figures  la),  b)  and 
c). 

TSM  is  created  by  rotating  about  the  X  axis  first,  then  about  Z,  then 
about  Y,  e.g., 

8)  TSM  =  Ry*RB*Rx 


1.  Matrices  TSM  and  T6S  are  named  to  conform  to  notations  used  in  Estes, 
1983. 
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An  older  version  of  the  FIT  program,  documented  in  Estes  (1983) ,  used  a 
different  order  of  rotation,  about  the  Z  axis  first,  about  the  new  I  axis 
second,  and  finally  about  the  new  Z  axis.  It  is  now  known  that  this  Z-X-Z 
rotation  order  fails  to  adequately  resolve  the  first  and  third  angles  when 
they  are  large,  and  so  the  present  order  of  rotation  was  instituted,  with 
successful  results. 

The  relation  of  euler  angles  6x>  £y»  cz  to  'roll,  pitch,  yaw*  notation 
is  dependent  on  the  spacecraft  axis  designations.  For  exaaple,  since  the 
MAGSAT  Z-axis  is  pointed  in  its  along- track  direction,  ez  is  roll, 

€j (radial)  is  yaw,  and  £x (cross  track)  is  pitch. 

The  FIT  program  calculates  some  field  quantities  in  earth-fixed 
cartesian  coordinates.  The  coordinate  origin  is  at  earth’s  center;  the  X 
axis  points  along  0*  longitude;  the  Y  axis  points  along  90°  meridian;  and 
the  Z  axis  points  along  the  geographic  north  pole.  Information  is  therefore 
required  on  the  relation  between  the  corrected  spacecraft  measurement 
vector,  Bspc,  and  its  analog  in  earth-fixed  coordinates  Bef,  for  every  data 
point.  This  information  is  contained  in  transformation  matrix  TGS: 

9)  Bef  =  TGS  *  Bspc 

The  TGS  matrix  itself  is  an  approximation  computed  from  the  formula: 


10)  TGS  =  f-nx  -rx  vx1 

I  I 

I  -Hy  -ry  Vy  I 

I  I 

l-nz  -rz  v^J 


where:  rx  =  cos^cosX,  ry  =  cos^sinX,  rz  =  sin^ 

nx  =  cos^nc°sXn,  ny  =  cos^QsinXQ,  nz  =  sin^n 

▼x  =  nyr*  -  ryna,  Vy  =  nzrx  -  rznx,  vz  =  nxry  -  rxny 


X,  and  XQ  are  defined  as: 

^  =  Geocentric  latitude 

X  =  East  longitude 

=  -8. 74* ,  the  inclination  of  the  vector  normal  to  the  orbit 

Xq  =  X  t  arccos[-tan(^n)*tan(^)] ,  where  +  is  used  for  a 

descending  orbit  (N  to  S)  and  -  for  an  ascending  orbit. 
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It  should  be  noted  that  this  equation  is  not  accurate  for  latitudes  greater 
than  about  75°.  DMSP  data  between  75°  and  81.26°  were  therefore  not 
utilised  in  subsequent  analyses.  Derivation  of  this  restriction  and  of  the 
TGS  transformation  eleaents  themselves  may  be  found  in  Appendix  fi. 

Combining  equations  #5  and  #0  yields: 

11)  Bef  =  TGS  *  TSM  *  TCAL*  (Bspu  -  bias) 

Bspu  and  TGS  are  read  or  computed  from  input  data,  and  TSM,  TCAL  and 
bias  are  solved  for  in  execution  of  the  FIT  program. 


IV.  Noise  Sources  in  the  Magnetic  Data. 


The  DMSP  data  were  examined  initially  by  Rich  (1984) ,  who  found  three 
sources  of  magnetic  noise.  The  first  two  are  high  frequency  sinusoidal 
signals  with  periods  of  0.576  and  3.456  seconds.  These  are  caused  by  the 
rotating  X-ray  scanner,  designated  the  SSB/S  instrument,  which  is  mounted  10 
to  15  inches  from  the  magnetometer  sensors  and  generates  a  small  magnetic 
field.  These  high  frequency  noise  sources  are  of  magnitude  less  than  about 
30  nT  and  are  not  a  concern  in  the  present  study. 

The  third  noise  source  found  by  Rich  is  the  operation  of  the  satellite 
torquing  coils.  These  are  turned  on  for  durations  of  about  4  minutes  at 
various  times  throughout  the  DMSP  mission.  When  the  coils  are  on,  the 
magnetic  field  data  is  offset  by  a  constant  level  shift  of  3000  to  14,000 
nT.  This  type  of  noise  is  screened  out  by  the  "gross  outlier"  criterion  for 
reducing  the  data  (see  section  VII) . 

A  foe  'th  known  noise  source  consists  of  fields  in  the  100  -  150  nT 
range  which  result  from  turning  on  transmitters  and  tape  recorders  when  over 
a  tracking  station.  It  is  assumed  that  most  data  so  affected  will  be 
eliminated  in  the  various  outlier  tests  described  in  section  VI. 

A  noise  source  not  discussed  by  Rich (1984)  was  discovered  by  examining 
orbital  plots  of  residual  data,  i.e.  data  which  have  had  a  preliminary  field 
model  subtracted.  These  residual  data  show  strong  periodic  trends  with 
amplitudes  of  up  to  70  nT,  after  other  corrections  were  applied.  An  in- 
depth  discussion  of  this  periodic  noise  and  its  removal  is  found  in  section 
VII. 
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V.  Observatory  and  Survey  Data  used  in  Modeling. 

Data  fros  land-based  magnetic  observatories  and  surveys  are 
instrumental  in  augmenting  the  satellite  data  for  modeling  purposes. 
Observatories  have  permanent  locations  and  are  useful  in  measuring  the  field 
behavior  over  an  extended  time  period.  They  usually  measure  declination, 
plus  the  horisontal  and  vertical  field  components,  from  which  X,  Y,  Z 
magnetic  components  (North,  East,  and  Down  in  a  geodetic  system)  may  be 
derived.  The  observatory  data  used  in  the  DMSP  studies  were  obtained  from 
the  World  Data  Center  in  Boulder,  CO.  They  are  annual  means  of  daily 
measurements,  so  each  falls  on  a  half  year  (e.g.,  1976.5).  Since 
observatory  biases  (see  next  section)  were  solved  for,  only  those  stations 
with  at  least  3  measured  components  and  3  years  of  data  were  retained. 

The  time  distribution  of  observatory  data  used  in  this  study  is: 


Table  2 


Time  period 

1979  -  1980 
1680  -  1981 

1982  -  1983 

1983  -  1984 

1984  -  1985 

1985  -  1986 

1986  -  1987 


#  of  stations 

143 

151 

155 

149 

140 

91 

1 


The  geographic  distribution  of  observatories  is  roughly  equal  to  that 
displayed  in  Langel  et.al.  (1988b). 

Survey  data  include  aeromagnetic  data  (Project  Magnet  and  a  Caribbean 
survey)  and  land  data.  The  Project  Magnet  data  are  3-component  vector 
recordings,  the  Caribbean  data  are  scalar  values  only,  and  the  land  data  may 
be  declination,  inclination,  horisontal,  north  (X),  east  (Y) ,  vertical  (Z) 
or  scalar  readings.  The  aeromagnetic  data  were  averaged  over  a  220  km 
distance  prior  to  the  least-squares  fitting  procedure,  in  order  to  reduce 
data  volume,  attenuate  crustal  anomalies,  and  obtain  an  estimate  of  data 
quality  for  data  weighting  in  the  fitting  procedure.  The  averaging  method 
is  described  in  detail  in  Langel  et.  al.  (1988b). 

The  land  data  were  checked  for  quality  by  comparing  to  a  preliminary 
field  model,  and  assigned  an  error  estimate  via  a  geographic  binning 
procedure  (see  Langel  et.  al.,  1988a  for  description  of  binning  process). 

The  data  have  extensive  geographic  coverage;  the  temporal  coverage  is 
described  in  Table  3. 
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Table  3. 


Magnet  data  Caribbean  Land  data 


Time 

X 

Y 

_ Z 

B 

D 

I 

H 

X 

Y 

Z 

B 

1980 

— 

— 

— 

1285 

110 

761 

... 

... 

521 

159 

1981 

306 

306 

306 

— 

105 

13 

88 

- — 

— 

71 

1982 

436 

436 

436 

— 

132 

39 

88 

— 

— 

68 

34 

1983 

278 

276 

276 

— 

105 

— 

124 

— 

— 

44 

60 

1984 

- — 

— 

— 

142 

92 

— 

71 

8 

7 

71 

3 

1985 

— 

— 

— 

— 

329 

89 

151 

— 

— 

105 

151 

Total 

1018 

1018 

1018 

142 

2048 

251 

1283 

6 

7 

880 

407 

VI.  Calibration  Procedure 


Fluxgate  magnetometers  are  subject  to  scale  and  zero-level  drifts  due 
to  changes  in  mechanical  properties  and  electronic  components.  For  MAGSAT, 
these  mere  solved  for  in  flight  by  comparison  with  the  on-board  scalar 
magnetometer  data.  Since  no  scalar  magnetometer  was  available  for  DMSP,  an 
attempt  was  made  to  use  data  from  magnetic  observatories  as  a  1 standard” 
against  which  to  calibrate  the  DMSP  data.  Calibration  is  accomplished  by 
solving  for  the  matrix  TCAL  (see  equation  #5)  in  the  modeling  procedure.  In 
this  process  the  'standard"  data  constrain  the  main  field  solution 
sufficiently  to  permit  the  TCAL  solution. ' 

Two  methods  of  calibration  were  developed,  both  of  which  use 
observatory  data.  The  use  of  such  data  in  a  main  field  model  necessitates 
taking  into  account  the  presence  of  crustal  fields,  fixed  measurement 
biases,  and  any  other  field  not  represented  by  the  model.  Langel  et.  al. 
(1082)  accomplished  this  by  designating  the  sum  of  such  fields  Bo,  and  the 
measured  field  at  the  observatory  B.  Bn  is  called  the  'bias”  at  a  given 
observatory .  B  equals  Bg  +  Bi,  where  Bi  is  the  field  from  the  spherical 
harmonic  model. 

Is  the  first  method  of  calibration,  Bm  is  determined  for  each 
observatory  prior  to  the  spherical  harmonic  solution,  under  the  assumption 
that  the  GSFC(12/83)  model,  epoch  1080,  is  highly  accurate  (Langel  and 
Bstes,  1985).  Bn  for  any  observatory  is  taken  to  be  the  measured  annual 
mean  nearest  1980  (generally  at  1980.5)  minus  the  field  value  computed  from 
the  GSFC (12/83)  model,  truncated  at  degree  10.  (For  consistency,  the 
truncation  level  is  the  same  as  that  used  for  models  derived  with  DMSP 
data.)  If  an  observatory  has  no  data  for  1980  then  the  bias  is  computed  at 
a  later  year  (1981.5  or  1982.5),  and  that  observatory  is  assigned  a  larger 
error  estimate  in  the  subsequent  spherical  harmonic  fit. 
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Once  Bb  lor  each  observatory  is  thus  obtained,  it  is  assumed  to  remain 
constant  through  tiae.  Bn  is  then  subtracted  from  measured  obserratory 
ralues  around  the  epoch  of  the  DMSP  data  being  calibrated.  The  corrected 
obserratory  data  are  then  interpolated  to  the  arerage  epoch  ol  the  DMSP  data 
and  a  field  model  is  generated.  The  DMSP  data  are  assigned  a  200nT  weight, 
compared  with  SnT  for  the  observatories .  In  the  spherical  harmonic 
solution,  satellite  biases,  Euler  angles,  and  the  elements  of  TCAL  are 
allowed  to  rary  in  addition  to  the  main  field  coefficients,  to  obtain  the 
calibration  constants. 

In  the  second  calibration  procedure,  Bn  is  solred  for  simultaneously 
with  the  main  field  coefficients.  This  solution  requires  that  data  from  at 
least  three  epochs  exist  for  each  obserratory  (Langel  et  al ,  1082) .  Thus 
all  obserratory  data  from  1080  through  the  DMSP  time  period  are  used,  as 
opposed  to  using  a  single  epoch  of  extrapolated  data,  as  in  the  first 
procedure.  The  observatories  and  DMSP  data  are  giren  weights  based  on 
residuals  to  prerious  fits  (generally  25  nT  for  observatories  and  40-80  nT 
for  DMSP  data) .  The  epoch  of  the  model  is  set  to  1080.  Both  main  field  and 
tiae  terms  are  computed,  along  with  the  DMSP  satellite  biases,  Euler  angles 
and  TCAL  elements.  The  GSPC (12/83)  information  matrix  is  used  as  a  priori 
information,  which  "anchors"  the  constant  solution  at  1080.  Observatory 
biases  are  solved  for,  i.e.  allowed  to  vary,  in  this  procedure,  versus  being 
held  fixed  in  the  first  procedure. 

This  second  procedure  gives  calibration  results  very  similar  to  the 
first  procedure.  However,  the  second  procedure  is  more  versatile,  since  it 
can  be  used  to  calibrate  the  DMSP  data  over  its  entire  time  span,  whereas 
the  first  procedure  can  only  be  used  to  calibrate  data  from  a  single  epoch, 
since  no  secular  variation  is  included  in  the  model.  Both  of  the  methods 
rely  on  the  GSFC  (12/83)  model  for  determining  the  observatory  biases,  with 
the  implicit  assumption  that  it  is  a  highly  accurate  model  and  that  the 
biases  remain  constant  over  the  tiae  span  1980  -  1985. 


VII.  Automated  Processing  Procedure  for  DMSP  Data. 


Preliminary  field  modeling 

A  test  model  was  generated  from  DMSP  data  at  epoch  1984.04.  The  giO 
term  from  this  model  equaled  -29,900.4  nT.  A  model  derived  from  observatory 
data  at  the  same  epoch  yielded  g]0  equal  to  -29,883.4.  The  closeness  of  the 
two  terms  suggests  the  apparent  adequacy  of  the  DMSP  data  for  main  field 
modeling.  A  calibration  of  the  1984.04  data,  using  the  second  procedure 
described  in  the  previous  section,  was  also  executed  with  the  following 
results:  SL1  *  0.9955,  SL2  =  0.9996,  SL3  =  1.0025.  The  nearness  of  these 
values  to  unity  again  suggests  that  the  DMSP  magnetometer  measured  the 
magnetic  field  accurately  for  that  selection  of  data.  These  results 
indicated  that  DMSP  data  might  be  useful  for  main  field  modeling,  in  spite 
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of  the  large  spacecraft  fields.  All  of  these  studies  were  conducted  using 
only  a  few  days  of  data.  On  the  basis  of  these  results,  it  was  decided  to 
proceed  with  a  larger  quantity  of  data. 

The  preliminary  field  aodel  was  removed  froa  January  14-18  DMSP  data  to 
create  residual  data.  Upon  examination,  these  orbits  of  data  showed  strong 
periodicities.  A  spectral  decomposition  of  the  data  revealed  noise  sources 
with  periods  equal  to  the  orbit  period  (100  minutes)  and  subharmonics  of 
1/2,  1/3  and  1/4  of  the  orbit  period.  Figure  2  displays  a  typical  residual 
orbit  froa  this  time  period.  The  X  and  Y  components  most  clearly 
demonstrate  this  periodic  noise.  Figure  3  is  the  associated  spectrum. 

Peaks  in  the  spectrum  display  these  dominant  noise  periods  quite  noticeably. 
As  will  be  seen  below,  one  cause  of  the  periodic  noise  is  the  need  for 
adjustment  of  the  Buler  angle  values  in  the  TSM  matrix  in  equation  5) .  The 
other  causes  of  this  periodic  noise  is  unknown.  The  peak-to-peak  amplitude 
of  the  noise  is  about  300  nT  before  Euler  angle  correction  and  about  50-70 
nT  after  that  correction. 


Automated  procedure  description 

A  four-step  procedure  was  followed  to  remove  data  spikes  and  periodic 
noise  froa  the  DMSP  data.  The  steps  are  as  follows:  1)  Fit  a  span  of  DMSP 
data,  covering  several  days,  with  a  preliminary  field  aodel.  Subtract  the 
field  model  to  get  residual  data.  Reject  data  points  above  75°  absolute 
latitude,  and  reject  'gross  outliers”,  i.e.,  residual  data  with  absolute 
values  greater  than  a  specified  cutoff.  Fit  the  residual  data  with  a  spline 
function  and  reject  points  which  deviate  more  than  2  standard  deviations 
from  that  function.  2)  Add  residual  data  which  is  not  rejected  back  to  the 
preliminary  field  model.  Then  fit  a  new  field  model  to  this  data  with  epoch 
equal  to  the  average  time  of  that  data  span.  Solve  for  constant  main  field 
coefficients,  magnetometer  angle  adjustments  and  biases.  3)  Correct  the 
original  data  with  the  angle  and  bias  solutions.  Use  the  computed  field  from 
step  #2  to  re-create  the  residual  data,  and  re-do  step  #1,  i.e.  reject  gross 
outliers  and  spline  outliers.  4)  Fit  a  Fourier  function,  which  is  composed 
of  the  4  dominant  noise  periods  (25  minutes,  33  minutes,  50  minutes  and  100 
minutes)  in  a  least-squares  manner  to  the  residual  data.  Reject  outliers 
according  to  the  Fourier  fit,  using  the  2a  criterion  as  for  the  spline  fit. 
Then  subtract  the  Fourier  function  from  the  data.  Add  the  result  back  to 
the  computed  field  model  from  step  #2,  to  create  the  final,  corrected  data 
set. 


Step  4)  is  somewhat  ad  hoc.  Such  periodic  variations  could  arise  from 
source  corrections  we  have  either  overlooked  or  been  unable  to  apply.  For 
example,  comparison  of  Figures  4a  and  4b  shows  that  much  of  the  large 
periodic  oscillation  results  from  unadjusted  Euler  angles.  It  both  is  more 
meaningful  and  reliable  to  correct  the  euler  angles  than  to  remove  the 
variations  via  the  Fourier  fit.  For  this  reason  the  ad  hoc  Fourier  fit 
correction  is  applied  last. 
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To  reiterate,  the  four  step  process  corrects  the  data  in  several  ways: 
a)  It  rejects  indiridual  bad  points  according  to  a  field  aodel,  b)  It 
rejects  points  above  78°  which  hare  an  inaccurate  attitude  determination,  c) 
It  corrects  the  data  with  magnetometer  rotations  and  biases  solved  for  in  a 
least-squares  fit,  d)  It  rejects  points  which  are  outliers  to  spline  and 
Fourier  function  fits,  and  e)  It  removes  periodic  noise  from  the  data  by 
subtracting  out  the  Fourier  function. 

Figure  4  shows  the  same  profile  from  Figure  2,  after  it  has  undergone 
the  data  cleaning  process.  Most  of  the  periodic  noise  is  gone,  and  the 
major  spikes  and  outliers  have  been  removed. 

Figure  5  is  a  world-wide  distribution  plot  of  3  days  of  processed  DMSP 
data.  Note  that  no  data  exists  polewaxdf  of  ±75°  latitude.  Also,  every 
orbit  contains  a  no-data  sone  approximately  15  degrees  long.  These  points 
have  been  rejected  according  to  the  'gross  outlier*  criterion,  and  evidently 
are  positions  in  the  orbit  where  the  torquing  coils  were  turned  on. 

A  diagnostic  of  the  quality  of  a  data  set  and  the  geomagnetic  model 
based  on  that  data  set  is  the  'spectrum'  from  that  field  model.  Consider 
the  quantity: 

13)  R„  =  (n+1)  E  [(gn®)2  +  (hn®)2] 

m=0 

Rn  is  the  mean-square  value  over  the  Earth's  surface  of  the  magnetic  field 
intensity  produced  by  harmonics  of  the  nth  degree.  Models  derived  from  the 
MA6SAT  data  are  considered  to  be  the  nearest  to  noise-free  of  all  available 
models.  The  presence  of  noise  at  any  degree  will  increase  the  value  of  Rn. 
Figure  6  shows  two  such  spectra  from  DMSP  models  compared  with  that  for  a 
MA6SAT  aodel  (Langel  and  Estes,  1982).  The  uncorrected  DMSP  spectrum,  from 
January  14-18  data,  begins  to  deviate  from  the  MA6SAT  spectrum  at  degree  9, 
whereas  the  spectrum  from  corrected  data  deviates  much  less,  and  not  until 
degree  10.  This  analysis  shows  that  the  removal  of  noise  sources  results  in 
an  improvement  of  the  character  of  the  output  field  models.  It  also  shows 
that  the  DMSP  data  are  considerably  noisier  than  the  MAGSAT  data,  even  when 
corrected. 


VIII.  Individual  Epoch  DMSP  Field  Models 


The  correction  procedure  was  applied  to  15  sub-sets  of  DMSP  data,  each 
containing  several  days  of  data.  Subset  epochs  ranged  from  January,  1984 
through  November,  1985.  Each  data  set  was  chosen  from  a  magnetically  quiet 
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period  as  determined  by  the  world-wide  Kp  index.  Results  of  step  2,  which 
soiree  for  the  field  model  and  magnetometer  adjustment  parameters,  are 
summarised  in  Table  4: 


Table  4. 


Date 

gl° 

gl1 

hi1 

Cx 

«y 

<« 

BIAS1 

BIAS2 

BIAS3 

yrs 

nT 

nT 

nT 

deg 

deg 

deg 

nT 

nT 

nT 

84.02 

-29895 

-1927 

5522 

-.038 

-.449 

0.005 

12.0 

2.6 

-1.7 

84.05 

-29893 

-1928 

5532 

-.065 

-.446 

0.006 

0.2 

-.5 

1.4 

84.21 

-29887 

-1935 

5523 

-.114 

-.459 

-.004 

7.8 

2.7 

-8.7 

84.34 

-29872 

-1925 

5516 

-.172 

-.457 

-.013 

11.4 

-2.0 

-12.5 

84.47 

-29860 

-1922 

5514 

-.138 

-.451 

-.016 

2.4 

-7.6 

-11.2 

84.63 

-29866 

-1932 

5503 

-.063 

-.474 

-.009 

-11.9 

-12.1 

-3.8 

84.71 

-29866 

-1927 

5505 

-.092 

-.478 

-.006 

-18.4 

-9.5 

-2.4 

85.05 

-29857 

-1918 

5496 

-.050 

-.496 

-.022 

-91.2 

-56.7 

20.0 

85.34 

-29856 

-1910 

5492 

-.129 

-.471 

-.013 

-91.2 

-67.0 

4.1 

85.45 

-29838 

-1920 

5494 

-.125 

-.467 

-.017 

-93.6 

-68.3 

1.3 

85.46 

-29843 

-1916 

5491 

-.130 

-.472 

-.020 

-86.7 

-65.9 

1.3 

85.60 

-29842 

-1915 

5495 

-.098 

-.475 

-.009 

-100.4 

-74.2 

4.6 

85.75 

-29847 

-1908 

5490 

-.074 

-.520 

0.013 

-112.6 

-72.3 

7.3 

85.82 

-29843 

-1914 

5484 

-.081 

-.511 

0.021 

-113.1 

-83.7 

12.8 

85.90 

-29832 

-1905 

5489 

-.030 

-.518 

0.032 

-110.0 

-57.5 

25.7 

Figures  7a)  through  7e)  are  derived  from  Table  4.  They  display 
solutions  for  g]0,  gj.1,  hil,  the  three  Euler  angles,  and  the  three  biases 
for  each  DMSP  data  set  throughout  time.  The  main  field  coefficients 
decrease  in  magnitude  with  time  as  expected  from  earlier  models,  but  the 
trend  is  not  smooth.  This  could  indicate  that  the  data  sub-sets  have 
marginal  geographic  distribution,  or  that  the  DIISP  data  are  not  sufficiently 
stable  over  time.  The  Euler  angle  solutions  are  fairly  consistent,  with  ey 
(yaw)  varying  slowly  from  -.44  to  -.52  degrees,  Cx  (pitch)  averaging  about 
-.1  degrees  and  (roll)  averaging  about  sero.  The  bias  values  show  a 
noticeable  break  between  September,  1984  and  January,  1985,  most  strongly  in 
X  and  Y.  Biases  at  January,  1985  depart  sharply  from  the  previous  bias 
trend  in  all  three  components.  This  jump  is  evident  in  the  biases  only,  and 
its  cause  is  uncertain.  One  possible  explanation  is  that  on  30  October, 
1984,  the  solar  array  panel  was  rotated  90°.  This  could  result  in  a  changed 
contribution  to  the  bias  field  from  the  solar  array  since  both  its  position 
and  its  total  current  were  changed.  Another,  though  less  likely,  spacecraft 
change  that  could  contribute  to  the  bias  change  is  that  on  7  November,  1984, 
the  skew  momentum  wheel  was  reset  so  that  it  drew  100  ma  less  current. 

The  bias  values  in  Table  4  are  part  of  the  value  of  the  vector 
parameter  biaa  to  be  used  in  equation  5),  i.e.  they  are  a  small  time 
dependent  correction  to  be  applied  in  addition  to  the  large  bias  values  of 
equation  4) .  A  small  further  correction  is  derived  in  section  IX. 


-14- 


H.  Field  Models  Using  All  DMSP  Data. 


Final  data  processing  steps 

Alter  the  DMSP  data  sub-sets  sere  corrected  according  to  the  euler 
angles  and  biases  in  Table  4,  they  were  prepared  lor  field  node ling  using 
the  entire  combined  set.  Before  concatenation,  data  froa  each  tiae  period 
were  first  exaained  according  to  3-hour ly  Kp  indices  and  hourly  DST  values. 
If  the  Kp  index  was  greater  than  2+  or  the  average  DST  value  was  greater 
than  lOnx  or  less  than  -30nT,  then  that  3  hour  section  of  data  was  rejected. 
Next,  data  froa  each  period  were  sorted  into  equal  area  geographic  bins  and 
randoaly  deleted  in  each  bin  until  the  best  approximation  to  the  following 
distribution  was  attained:  3  vector  (i.e.,  3  components  each)  points  per 
bin  for  dip-latitudes  less  than  30  degrees,  and  9  scalar  points  per  bin  for 
dip-latitudes  greater  than  30  degrees.  This  distribution  is  based  on  work 
done  by  Lowes  and  Martin  (1987) ,  which  indicates  that  vector  data  are 
required  at  low  magnetic  latitudes  to  avoid  the  'perpendicular  error 
effect”,  while  scalar  data  at  high  latitudes  are  more  effective  than  vector 
measurements,  because  they  are  more  accurate,  especially  in  the  presence  of 
field-aligned  currents. 

After  binning,  the  15  different  time  periods  were  merged  together  into 
a  single  file  and  entered  into  a  format  suitable  for  input  into  the  FIT 
program.  The  final  data  set  consisted  of  7750  X  and  Y  measurements,  7768  Z 
data,  and  16097  B  data. 

All  of  the  DMSP  data  were  then  calibrated  according  to  the  second 
procedure  in  section  VI.  The  field  model  for  this  calibration  was  of  degree 
eleven  in  the  constant  field,  first  derivative  (secular  variation  or  SV) 
terms,  and  second  derivative  (secular  acceleration  or  SA)  terms.  The 
terminology  used  herein  to  designate  such  a  model  is  (N1,N2,N3),  where  N1  is 
the  degree  of  expansion  of  constant  spherical  harmonic  terms,  N2  the  degree 
of  expansion  of  first  derivative  terms  and  N3  the  degree  of  second 
derivative  terms.  The  calibration  model  is  thus  described  as  (11,11,11). 
Note  that  N1^N2>N3.  The  GSFC (12/83)  model  was  used  as  a  priori  information, 
the  epoch  was  1980,  and  the  calibration  solution  included  satellite 
biases,  Euler  angles,  and  TCAL  elements.  In  the  solution,  the  Euler  angles 
were  negligible,  the  bias  vector  equaled  (4.8nT,  -0.51nT,  0.65nT),  and  the 
slope  solution  equaled:  SL1  =  0.99994,  SL2  =  0.99960,  SL3  =  1.00119. 

These  calibration  results  were  applied  to  the  data  before  creation  of 
the  final  field  models,  and  these  are  the  final  corrections  to  be  applied  in 
accord  with  equation  5).  Thus,  in  equation  5),  the  value  of  bias  is  the  sum 
of  this  bias  value  and  those  found  in  Table  4;  the  TCAL  matrix  is  derived 
from  these  values  of  SL1,  SL2,  and  SL3;  and  the  TSM  matrix  is  derived  from 
the  Euler  angles  of  Table  4  with  no  additional  correction. 
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Use  of  DB-2  data  in  models 


The  DMSP  data  were  combined  in  subsequent  fits  with  3-coaponent  and 
scalar  magnetic  field  data  iron  the  Dynamics  Explorer  2  (DB-2)  spacecraft. 

A  complete  description  of  these  data  is  given  in  Ridgway  (1988) ,  but  is 
summarised  here.  19,600  rector  DE-2  magnetic  observations,  between  10/81 
and  1/83,  were  made  available  to  the  Geology  and  Geomagnetism  branch, 

Goddard  Space  Flight  Center.  Since  the  primary  mission  of  DE-2  was  to  map 
field-aligned  currents  at  high  latitudes,  the  data  were  concentrated  in 
polar  regions.  The  data  were  also  not  well  distributed  temporally,  being 
concentrated  at  epoch  1982. 

These  problems  necessitated  implementation  of  a  binning  procedure  on 
the  data  and  decimation  within  the  bins,  to  distribute  the  data  more  equally 
both  temporally  and  geographically.  This  procedure  cut  the  number  of  data 
to  10,600  vector  observations,  improved  the  temporal  distribution,  reduced 
data  concentration  at  the  poles,  and  improved  data  quality  by  rejecting  bad 
points . 

The  DE-2  magnetometer  was  boom-mounted,  in  order  to  minimize  spacecraft 
fields.  Accuracy  of  the  spacecraft  attitude  determination,  however,  was 
only  0.5°-0.7°  (Langel  et.  al.,  1988),  resulting  in  bias  and  rotational  off¬ 
sets  in  the  three  magnetic  components.  Attitude  corrections  to  the  data 
were  empirically  computed  using  the  FIT  software  (see  section  III) ,  which 
significantly  reduced  residuals  from  calculated  field  models  when  applied. 
The  DE-2  data  were  also  calibrated  to  ground  observatories  in  the  same 
manner  as  the  DMSP  data.  Results  from  the  attitude  correction  and  the  cali¬ 
bration  are  as  follows: 

Euler  angles  :  ex  =  0.220°,  ey  =  -0.091°,  ez  =  0.411° 

Biases  :  biasi  =  -10.8nT,  bias?  =  6.9nT,  biasg  =  2.8nT 

Calibration  slopes  :  SLi  =  0.99928,  SL2  =  0.99976,  SL3  =  1.00011 

(For  biases  and  slopes,  along- track  =  1,  radial  =  2,  and  cross-track  =  3) 

Application  of  these  corrections  improved  the  data  fit  to  all 
subsequent  geomagnetic  field  models.  The  data  residuals,  relative  to  the 
DE-2  field  model  at  1982  (excluding  external  field  solution),  are: 

<7X  =  44  nT 
ffy  =  55  nT 
=  48  nT 
a\)  =  24  nT 

These  results  show  that  the  scalar  magnitude  data  (residual  <7b)  is  fit  much 
better  than  the  vector  component  data.  Later  models,  whether  using  DE-2 
data  alone  (see  Ridgway,  1988)  or  DE-2  plus  DMSP  data  (see  following 
section)  were  generally  most  successful  when  utilizing  only  the  scalar  DE-2 
data. 
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Initial  Models 


The  first  nodal  using  all  of  the  DMSP  data,  called  the  MQDMSDE  model, 
was  of  degree  (11,11,11),  neaning  main  field  static  terms,  linear  time 
terms,  and  2&d  order  time  terms  up  to  degree  11.  The  *MR  of  MQDMSDE  stands 
for  MAGSAT  a  priori  information,  from  the  GSFC(12/83)  model.  The  "0*  is 
for  observatory  annual  means  extending  from  1979.5  through  1985.5.  The  "DM* 
is  for  DMSP  data,  with  rector  data  at  low  latitudes  and  scalar  magnitudes 
abore  30°  dip  latitude.  The  *S*  stands  for  world-wide  surrey  data  of  all 
measurement  types,  which  were  weighted  with  a  correlated  weight  matrix  (see 
Langel  et.  al.,  1989).  The  *DB*  stands  for  DE-2  data,  dirided  into  rector 
and  scalar  data  just  as  for  the  DMSP  set.  The  model  was  created  in  steps,  in 
the  order  of  the  letters.  Also  solred  for  were  obserratory  biases,  three 
external  field  coefficients  of  the  first  order,  and  the  4  DST  multipliers 
(one  for  glP,  q]0,  qj.1  and  s]!)  .  The  model  epoch  was  1980,  and  at  1980  the 
model  was  hearily  constrained  by  the  MAGSAT  data,  with  the  other  data 
primarily  determining  the  temporal  rariation  terms. 


Examination  of  secular  rariation  terms 


The  secular  variation  results,  for  MO,  MODM,  MODMS,  and  MODMSDE,  all  of 
degree  (11,11,11)  at  epochs  1980,  1985,  and  1990,  are  displayed  in  Table  5: 

Table  5. 


Model 

Epoch 

*i° 

ftl1 

El1 

&2° 

S21 

R21 

1980 

25.7 

8.0 

-28.0 

-19.8 

5.2 

-3.3 

MO 

1985 

21.3 

9.2 

-20.4 

-17.0 

5.4 

-18.4 

1990 

17.0 

12.4 

-12.7 

-14.1 

5.6 

-33.5 

1980 

25.2 

6.2 

-20.48 

-17.2 

3.9 

-9.9 

MODM 

1985 

28.6 

8.9 

-20.43 

-18.7 

2.9 

-17.3 

1990 

32.1 

11.7 

-20.38 

-20.2 

1.9 

-24.6 

1980 

24.4 

9.5 

-23.3 

-16.8 

4.6 

-8.7 

MODMS 

1985 

29.0 

6.3 

-18.0 

-19.0 

2.0 

-18.0 

1990 

33.7 

3.1 

-12.8 

-21.2 

-0.5 

-27.4 

1980 

25.1 

15.2 

-19.5 

-15.7 

3.7 

-10.1 

MODMSDE 

1985 

28.5 

1.6 

-21.3 

-19.9 

2.8 

-16.9 

1990 

31.9 

-11.9 

-23.2 

-24.2 

2.0 

-23.8 
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The  behavior  of  these  models  is  not  consistent.  For  the  MO  model, 
decreases  with  time,  in  contrast  to  MODM,  MODMS  and  UODMSDE ,  where  it 
increases  fron  1980  to  1990.  For  the  MO  and  MODM  models  are  in  good 
agreement,  the  MODMS  is  substantially  different,  and  the  MODMSDB  shows  a 
totally  opposite  trend.  For  All,  both  MO  and  MODMS  increase  over  time,  at 
roughly  the  same  rate.  In  contrast,  the  MODM  All  remains  constant  through 
time,  and  the  MODMSDE  decreases. 

It  was  now  apparent  that  the  models  had  problems.  Besides  the 
inconsistent  secular  variation  results,  it  was  found  that  the  calculated 
correlations  between  the  SY  and  SA  coefficients  were  highly  negative, 
approaching  -1.  This  means  that  the  solution  has  strong,  undesirable 
interdependences  which  might  also  contribute  to  the  discrepant  secular 
variation  results.  The  parameter  space  was  selectively  scaled  back  to 
(11,10,5),  (11,8,5),  and  (11,8,0)  models.  Whenever  SA  terms  were  included  in 
the  solution,  high  negative  correlations  with  the  SY  terms  invariably 
resulted.  When  the  models  were  of  degree  (11,8,0),  most  of  the  high 
correlations  disappeared,  indicating  that  this  is  near  to  the  maximum 
parameter  space  which  can  be  supported  by  the  available  data. 

The  MO  model  was  also  examined.  It  also  showed  large  correlations, 
although  not  as  high  (±0.75),  and  more  scattered  in  the  solution  space. 

These  correlations  persisted  even  when  the  MO  solution  was  cut  to  a  (11,8,0) 
model,  indicating  that  available  observatory  data  alone  are  also  not 
sufficient  to  solve  for  a  model  with  this  size  parameter  space. 

Examination  was  also  made  of  the  secular  variation  results  from  the 
M0(ll,8,0),  (11,11,0)  and  MODM(11,10,5)  models.  &i°,  ftil  and  A^  terms  for 
these  models  are  plotted  in  Figure  8,  along  with  the  same  coefficients 
derived  from  observatory  first-difference  data  from  1940  through  1979.  If 
the  new  models  are  valid,  they  should  be  consistent  with  the  trend  of  the 
first  difference  plots.  The  &0  plots  generally  match  the  first-difference 
trend  for  the  MO  models.  The  MODM  model  is  somewhat  off  of  the  trend.  For 
£ll,  both  MO  (11,10,5)  and  MO  (11,8,0)  fall  roughly  in  the  expected 
coefficient  range.  MO  (11,11,0)  is  somewhat  out  of  the  range,  and  the  MODM 
model  is  definitely  off  of  the  trend.  For  the  All  plot,  all  of  the  MO 
models  and  also  the  MODM  model  are  consistent  with  the  trend  of  the  first- 
dif  f erence  coef  f icients . 

Clearly  the  different  ftll  results  in  Figure  8b  are  inconsistent. 
Agreement  of  the  M0(ll,8,0)  model  with  the  first-difference  trend  argues,  as 
above,  that  the  observatory  data  distribution  is  inadequate.  Addition  of 
the  DMSP  data  should,  in  principal,  greatly  improve  the  data  distribution 
and  permit  solution  for  more  parameters.  But  clearly  the  MODM  model  is  in 
serious  disagreement  with  the  first-difference  trend.  Adding  to  the 
confusion  is  the  apparent  agreement  of  the  MODM(11,10,5)  model  with  the 
MO (11, 11)  model. 
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To  this  point  DMSP  rector  data  has  been  utilized.  Since  rector  data  is 
■ore  error-prone  than  scalar  data,  it  was  decided  to  examine  the  possibility 
that  at  least  soae  of  the  model  problems  stemmed  from  systematic  errors  in 
the  DMSP  rector  data. 


M08DMDB  (11,8,5)  and  (11.11.51  scalar  models 

Errors  in  magnetometer  attitude,  either  random  or  systematic,  should 
not  affect  the  measured  scalar  magnitude.  Accordingly,  in  order  to 
eliminate  any  effects  of  strictly  attitude  errors,  in  the  next  series  of 
models  only  scalar  DMSP  and  DE-2  data  were  used,  in  addition  to 
obserratories  and  surrey  data.  The  data  included  13,121  DMSP  data  points  at 
an  arerage  epoch  of  1985,  and  5100  DE-2  points  with  an  arerage  epoch  equal 
to  1982.3.  The  satellite  data  were  added  to  a  base  model  consisting  of 
MAGSAT  a  priori,  obserratory  data,  and  surrey  data,  called  MOS.  The  epoch 
of  both  the  (11,8,5)  and  (11,11,5)  models  is  1980. 

In  both  resulting  models  the  secular  acceleration  terms  again  hare  high 
correlations  with  the  corresponding  secular  rariation  terms,  gi^,  gi*  and 
All  coefficients  for  each  model  are  displayed  in  Figure  9,  along  with 
coefficients  derired  from  the  first-differences.  (For  this  plot,  first- 
difference  coefficients  were  also  generated  from  1980-1985,  to  add  to  the 
1940-1979  data  used  in  Figure  8.)  Plotted  in  Figure  9  are  the  MO  (11,8,5), 
MOS  (11,11,5),  and  the  MOSDM  and  MOSDMDE  (11,8,5)  and  (11,11,5)  models,  at 
epochs  1980  and  1985. 

Examination  of  Figure  9  again  focuses  attention  on  gil.  The  gi°  and 
All  models  are  all  in  reasonable  agreement.  Models  of  gil  also  show  some 
sense  of  internal  agreement,  but  now  with  a  totally  unexpected  trend.  The 
ftl1  coefficient  decreases  sharply  but  steadily  from  1982  to  1985,  both  in 
the  first-difference  coefficients  and  in  the  models  utilizing  satellite 
data.  Models  with  observatory  and  survey  data  [MO (11, 8, 5)  and 
MOS (11, 11, 5) ] ,  which  have  SA  terms,  show  the  same  trend,  though  not  as 
pronounced.  This  is  highly  unexpected  behavior,  unless  some  sort  of 
geomagnetic  jerk  is  in  progress. 


MOSDMDE  (11.11, 5. 3)  models 

The  possibility  of  a  jerk  in  the  field  at  1982  prompted  a  test  of  a 
(11,11,5,3)  model,  which  provides  for  non-linear  fluctuation  of  the  secular 
variation  terms.  Only  scalar  DMSP  and  DE-2  data  were  used,  as  in  the 
previous  models.  This  model  was  produced  despite  the  high  correlations 
inherent  in  solving  for  an  overabundance  of  time  terms.  Once  MOSDM 
(11,11,5,3)  and  MOSDMDE  (11,11,5,3)  models  were  completed,  parameters  were 
also  cut  back  to  obtain  a  MOSDMDE  (11,9,5,3)  model.  This  proved  to  be  very 
similar  to  the  (11,11,5,3)  models. 
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Plots  of  the  secular  variation  coefficients  from  these  models  present 
no  consistent  picture.  For  and  the  models  with  satellite  data  are 

in  disagreement  with  the  first-difference  trends.  For  fill,  82®,  82 1  and 
the  trends  are  similar,  but  not  exact.  The  models  for  gll  snow  a  decreasing 
tend  after  1980,  but  the  first-difference  trend  differs  significantly  from 
the  trend  of  models  with  satellite  data.  Examination  of  first-difference 
plots  of  observatory  data  showed  no  conclusive  evidence  for  a  "jerk*  in  this 
period.  However  this  is  not  a  conclusive  result  since  the  time  span  of 
available  data  from  most  observatories  was  too  short  to  definitively 
determine  the  trend.  Our  conclusion  is  that  the  data  available  are  not 
adequate  to  discern  the  existence  of  a  jerk. 


Trends  of  the  spherical  harmonic  coefficients 

We  are  left  with  the  pussle  of  a  possibly  anomalous  behavior  of  the  &.1 
coefficient,  particularly  when  DMSP  data  are  combined  with  observatory  and 
surface  data.  It  is  not  yet  clear  that  the  DMSP  data  are  consistent  with 
the  surface  data,  or  that  the  DMSP  data  are  reliable  for  geomagnetic  field 
modeling. 

Figure  10,  which  plots  the  g]0,  g].l,  and  hil  coefficients,  not  the 
derivatives,  derived  from  the  IGHF  (1970,1975,1980),  the  DE-2  calibration 
(1982),  and  the  individual  DMSP  runs  from  Table  4,  shows  more  clearly  the 
nature  of  the  problem  with  the  models  based  on  the  DMSP  data.  For  the  g]0 
plot  (Figure  10A) ,  the  DMSP  results  fall  fairly  well  along  the  linear  trend 
of  the  IGRF  models.  The  DE-2  calibration  term,  at  1982,  is  slightly  high, 
but  not  significantly  off -trend.  The  gll  plot  (Figure  10B)  reveals  much 
more  about  the  DMSP  data.  All  of  the  gll  coefficients  derived  from 
individual  DMSP  fits  are  significantly  off  of  the  trend  of  the  IGRF  models, 
having  a  lower  magnitude  than  expected.  The  DE-2  gll,  however,  is  slightly 
higher  than  expected  indicating  a  discrepancy  between  the  two  data  sets. 

The  hil  term  (Figure  10C)  apparently  varies  in  a  non-linear  manner,  with 
both  the  DE-2  and  the  DMSP-derived  coefficients  falling  along  the  general 
trend  of  the  IGRF  models. 

Figure  11  allows  further  examination  of  the  behavior  of  the  computed 
gll.  Figure  11A  shows  the  individual  DMSP  fits,  along  with  the  DE-2  and 
DMSP  calibration  models.  The  DMSP  calibration  matches  the  individual  DMSP 
fits  very  well.  It  diverges  noticeably  from  the  DE-2  calibration.  Figure 
11B  displays  the  individual  DMSP  fits  along  with  the  (11,9,5,3)  MOSDMDE 
model,  the  (11,9)  MOSDMDE  model,  and  the  MOS  (11,9)  model.  This  plot  shows 
that  the  DMSP  models  depart  noticeably  from  the  models  based  only  on 
observatory  and  survey  data  (MOS  model) .  The  DMSP  data  require  a  lower  gll 
term.  When  combining  DMSP  data  with  observatory  and  repeat  data,  as  in  the 
M0SDMDE(11,9,5,3)  model,  the  gll  term  must  change  dramatically  with  time  to 
accommodate  all  data  types. 
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X.  Conclusions 


The  DMSP  body-mounted  magnetometer  data,  recorded  between  1984  and 
1988,  hare  been  examined  for  their  utility  in  geomagnetic  field  modeling. 

The  data  contain  major  sources  of  error,  as  indicated  by  large  data  spikes 
of  up  to  several  thousand  nT,  periodic  noise  of  several  hundred  nT,  and 
rotational  and  bias  errors  in  the  magnetometer  of  up  to  0.5  degree  and  one 
hundred  nT. 

The  data  spikes  are  corrected  by  a  gross-outlier  cutoff  relative  to  a 
preliminary  field  model  and  a  spline  fit.  Periodic  noise  is  removed  by 
fitting  the  data  residuals  with  a  4-component  Fourier  series.  Magnetometer 
attitude  errors  are  corrected  by  solving  for  rotational  offsets  using  the 
FIT  software,  and  the  biases  are  similarly  corrected.  The  total  magnitude 
of  the  data  is  adjusted  by  calibrating  it  against  ground  observatories. 

These  corrections  improve  the  data  quality,  as  indicated  by  residual 
plots  and  by  power  spectra  of  resultant  field  models.  However,  although  the 
field  models  show  magnitudes  roughly  within  the  realm  of  believability, 
closer  inspection  reveals  basic  flaws.  Any  models  using  Taylor  series  time 
representations  give  high  correlations  between  the  secular  variation  and 
secular  acceleration  parameters.  This  is  true  for  both  vector  and  field 
magnitude  data.  More  importantly,  the  models  derived  from  fits  to  the 
individual  data  sets  show  significant  differences  to  the  trend  of  previous 
IGBF  models  and  models  derived  from  DE-2  data.  These  differences  are 
manifested  particularly  in  the  gj.1  spherical  harmonic  coefficient.  In 
hindsight,  if  the  plots  displayed  in  Figure  11  had  been  examined  early  in 
the  analysis  process,  much  of  the  confusion  over  the  secular  variation 
results  could  have  been  avoided. 

Unless  further  methods  are  found  to  correct  these  discrepancies,  it  is 
concluded  that  the  data  at  present  are  not  suitable  for  main  field  modeling. 
There  are  several  possible  reasons  for  this.  First,  the  spacecraft  fields 
may  simply  be  so  strong  or  variable  that  any  correction  procedures  are 
destined  to  fail.  Second,  in  addition  to  hard  magnetic  sources  and  fixed 
amplitude  variable  magnetic  sources  (such  as  currents) ,  there  may  be  a 
substantial  amount  of  soft  magnetization  on  the  spacecraft.  This  would 
result  in  induced  fields  which  would  vary  with  time  and  position,  and  for 
which  we  do  not  at  present  know  how  to  correct.  Third,  the  timing  and 
position  error  on  the  data  furnished  to  Goddard  Space  Flight  Center  may 
result  in  systematic  errors.  Fourth,  the  automated  procedure  for  removing 
noise  from  the  data,  particularly  the  removal  of  Fourier  components,  may  be 
inadequate  and  could  possibly  introduce  error  into  the  resulting  data. 

Nothing  can  be  done  about  the  first  and  second  items  in  the  present 
data.  If  further  effort  were  to  be  expended  upon  the  DMSP  F-7  data,  it 
should  deal  with  items  three  and  four.  A  possible  approach  would  be  as 
follow:  a)  obtain  new  data  from  the  AFGL,  for  the  same  quiet  day  periods, 
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with  no  round-oil  in  tine  and  position  information,  b)  Eliminate  spurious 
data  by  plotting  all  of  the  data,  in  short  sections  (1  or  2  orbits) ,  and 
hand-selecting  data  segments  which  are  relatively  free  from  variable 
spacecraft  fields  due  to  equipment  being  turned  on  or  off.  Hopefully  this 
procedure  would  make  the  Fourier  correction  unnecessary. 

It  is  not  certain  that  these  procedures  would  be  able  to  salvage  the 
DMSP  F-7  data  for  main  field  modeling  purposes.  However,  the  exercise  of 
processing  the  data  has  been  instructive,  despite  difficulties  encountered, 
and  good  possibilities  exist  for  use  of  future  data. 

If  the  measurement  accuracy  can  be  sufficiently  upgraded,  future  DMSP 
data  could  make  a  significant  contribution  to  main  field  geomagnetism.  This 
is  because  the  present  need  is  for  long-term  operational  observatory-class 
measurements.  Such  data  would  not  only  provide  operational  tools  such  as 
field  models,  navigation  charts  and  submarine  magnetic  detection  baselines, 
but  would  provide  a  research  tool  desperately  needed  to  probe  the  nature  of 
the  geodynamo,  the  conductivity  of  the  mantle,  and  the  core/mantle  momentum 
exchange  mechanism.  Keys  to  upgrading  the  data  quality  are:  1)  placing  the 
magnetometer  on  the  end  of  a  boom,  2)  upgrading  the  attitude  determination 
at  the  magnetometer  location,  and  3)  increasing  the  accuracy  of  the 
magnetometer . 

It  is  our  understanding  that  step  3)  has  already  been  taken,  at  least 
in  part.  A  boom-mounted  magnetometer  would  significantly  reduce  periodic 

noise  and  bias  errors  in  the  data.  Combined  with  more  accurate  attitude 
determination,  such  an  addition  to  the  DMSP  magnetic  measuring  capabilities 
would  greatly  enhance  the  usefulness  DMSP  data  in  deriving  accurate  main 
field  models. 
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XI .  Appendix  A:  DMSP  Data  Format 


I.  Tape  parameters 

Tape  density:  6250  BPI 

Logical  record  length:  75  bytes 
Block  sise:  1875. 

Fixed  block,  ASCII. 


H.  Data  format 

The  data  are  arranged  as  one  ephemeris  record  syery  minute, 
followed  by  80  magnetometer  records  (one  per  second) . 

A.  Bphemeris  record 


location 

format 

description 

1-  4 

14 

Year. 

5-  8 

14 

Day  of  year 

9-14 

16 

Time  of  ephemeris,  seconds  U.T. 

15-18 

14 

Altitude  (nautical  miles) . 

19-28 

F10.2 

Geographic  latitude. 

29-38 

F10.2 

Geographic  longitude. 

?"  48 

F10.2 

Corrected  geomagnetic  latitude. 

49- w 

F10.2 

Corrected  geomagnetic  longitude. 

59-68 

F10.2 

Corrected  geomagnetic  local  time 

69-71 

3x 

Blank. 

72-75 

14 

Number  of  data  records  following 

ephemeris  (usually  =  60) . 


B .  Data  record 


location 

format 

description 

1-  8 

16 

Time  of  data  record  (seconds  UT) 

7-11 

5x 

Blank. 

12-29 

318 

X,  Y,  Z  magnetometer  counts  for 
first  of  20  samples  per  second. 

30-33 

4x 

Blank. 

34-51 

318 

X,  Y,  Z  magnetometer  counts  for 
eleventh  of  20  samples  per 
second. 

52-55 

4x 

Blank. 

58-75 

1012 

Ten  data  quality  flags. 
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TTT.  Appendix  B:  Derivation  of  DMSP  Transformation  Matrix,  TGS. 


For  every  DMSP  magnetic  reading,  known  position  information 
includes  the  geographic  latitude,  longitude  and  altitude  of  the 
satellite,  and  the  inclination  of  the  satellite  orbit  plane.  From  this 
information,  an  operator  may  be  derived  which  transforms  measurements  in 
the  local  satellite  coordinate  system  to  an  earth-fixed  cartesian 
coordinate  system  (Originally  derived  by  D.  Chinn  of  Science 
Applications  Research  Inc.). 

Figure  B1  displays  a  DMSP  satellite  orbit  superimposed  upon 
rotating  Earth-fixed  cartesian  coordinates.  The  inclination  of  the 
satellite  orbit  plane  with  the  earth’s  equatorial  plane  is  98.74 
degrees.  A  local  orthogonal  coordinate  system  may  be  assigned  to  the 
satellite,  with  one  axis  pointing  along  the  satellite  track,  another 
normal  to  the  satellite  track,  and  a  radial  axis  pointing  outwards  from 
the  earth-fixed  origin.  The  corresponding  unit  vectors  are  v,  n,  and  r, 
respectively . 

The  unit  radial  vector  r,  in  cartesian  coordinates,  equals: 


Bl) 

r  *  rxi  +  ryj 

where : 

B2a) 

rx  =  cos^cosX 

B2b) 

ry  as  cosjlsinX 

B2c) 

r*  =  sinjj 

and  where  ^  equals  the  latitude,  X  the  longitude  of  the  satellite. 

The  unit  normal  vector  n  may  be  described  in  an  analogous  ma.iner, 
except  that  In  replaces  f  and  XQ  replaces  X  in  equations  B1  and  B2.  ftn 
equals  a  constant  -8.74s  anywhere  on  the  orbit,  while  XQ  varies 
throughout  the  orbit.  XQ  is  solved  for,  given  ^  and  X,  by  utilizing  the 
fact  that  r  *  a  equals  zero: 

B3)  cos^sinXcos^nsioXn  ♦  cos^cosXcos^ncosXn  +  sin^sin^n  =  0* 

Dividing  by  cos^  and  cos^n  and  switching  tangents  to  the  other  side 
of  the  equation  yields: 

B4)  -tan^tan^n  =  cosXcosXn  +  sinXsinXQ. 

Utilizing  the  addition  formula  for  cosines  yields: 

B6)  -tan^tan^n  =  cos (X  -  Xn) 

or: 
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B6)  Xn  *  X  ?  arccos (  -tan^tan^n  )  (-  for  ascending  orbit, 

+  for  descending  orbit) . 

Once  Xn  is  solved  for,  then  r  and  a  are  fully  determined  by  B2.  ▼ 

is  derived  froa  the  relation: 

B7)  v  =  a  x  r 


The  transformation 
coordinates  is  thus: 

from 

cartesian 

(i,J,k)  to  local  (r,n,v) 

B8)  f  r  1  (tx 

1  1  1 

Ini  =  lnx 

1  1  1 

l  v  J  lvx 

r7 

\ 

nBl  * 

1 

(  i  1 
|  | 

n7 

1  J  1 

1  I 

y7 

1  1 
l  k  J 

The  MAGSAT-compatible  satellite  coordinates  (x-cross  track,  y- 
radially  down,  a -along  track)  correspond  to  -a,  -r,  and  v,  respectively, 
so  that  the  rotation  froa  cartesian  to  spacecraft  coordinates  is  thus: 


BO) 


(cross') 

I  I 

I  down  I 

I  I 

lalongj 


r-nx 

1 

y 

—nz1 

| 

rn 

i  i 

1 

l-rx 

-ry 

* 

-rzl  * 

131 

1 

1 

1  1 

'■▼x 

y7 

▼a  J 

IkJ 

The  TGS  matrix  is  a  rotation  from  spacecraft  to  cartesian 
coordinates,  and  so  equals  the  inverse  (=  transpose  for  a  rotation 
matrix)  of  the  above  equation: 


fll 

1  i 

r-nx 

I 

~rX 

1 

f  cross'! 

1  l 

1  1 

HI 

I  | 

1 

■  l-ny 

1 

-ry 

Ty! 

1  1 

Idown  1 

I  | 

UcJ 

1 

l-nz 

-rz 

VZJ 

1  1 

lalongj 

- TGS - 


Accuracy  of  derivation 


The  accuracy  of  the  TGS  derivation  depends  strongly  on  the  accuracy 
of  the  XQ  value  in  equation  B8) .  A  small  error  in  longitude  (X) 
determination  will  propagate  to  an  equal  error  in  XQ,  but  an  error  in 
latitude  (^)  will  propagate  to  a  larger  error  in  Xn.  This  relation  is 
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seen  by  taking  the  derivative  of  Xn  with  respect  to 

Bll)  dXn/d^  *  -K/(  cosf/cos2(fl  -  K2sin2(f)  )  with  K  =  tan (-8. 74) 
The  relation  of  dXn/d^  to  ^  is  as  follows: 


TAHTiR  B1 


*  dXn/d* 

81.28°  5 

80*  10.4 

78*  5.2 

78*  3.3 

75*  2.8 

If  l  is  known  to  within  0.2° ,  for  instance,  then  the  uncertainty  in 
Xn  would  equal  0.58,  0.87,  1.0  and  2.1  at  75°,  78°,  78°,  and  80°, 
respectively.  Because  of  this,  it  was  decided  not  to  accept  data  above 
75°  latitude. 

Alternative  formulation 


An  alternative  approach  to  defining  the  TGS  matrix  could  be  taken 
if  certain  DMSP  orbital  parameters  were  available.  The  following 
derivation  is  included  for  clarity  of  the  meaning  of  TGS.  Consider  a 
typical  satellite  orbit  in  fixed  inertial  space,  e.g.,  Figure  B2.  The 
X-Y-Z  coordinate  system  is  fixed  in  inertial  space  with  Z  along  the 
Barth’s  rotation  axis  and  X  along  the  vernal  equinox.  The  vernal 
equinox  is  the  ascending  node  of  the  Earth’s  orbit  around  the  sun.  The 
epoch  of  X  and  Z  must  be  specified.  Y  is  chosen  to  make  a  right  handed 
system.  For  our  purposes,  consider  the  plane  of  the  orbit  to  be  fixed 
in  the  X-Y-Z  system.  Define  the  spacecraft  coordinate  system,  X’-Y’-Z’ 
as  follows:  X’  will  be  positive  in  the  radial  direction,  Z’  will  be 
normal  to  the  orbit  plane,  with  Y’  chosen  to  make  a  right  handed  system. 
For  circular  orbits,  Y’  is  along  track.  The  angle  of  the  orbital  plane 
with  respect  to  the  equatorial  plane  is  called  the  inclination,  i;  the 
intersection  of  the  orbital  plane  and  the  equatorial  plane,  at  the  point 
where  the  satellite  is  going  north  is  the  ascending  node.  The  angle 
between  X  and  the  ascending  node,  Q,  is  called  the  right  ascension  of 
the  ascending  node.  The  satellite  position  in  the  orbit  is  then 
specified  by  its  angle  in  the  orbit  plane,  u,  measured  from  the 
ascending  node. 

To  transform  a  vector  or  position  from  the  X’  system  to  the  X 
system  requires  a  set  of  three  rotations:  (1)  a  rotation,  Ri,  around  the 
Z’  axis  through  the  angle  w;  (2)  a  rotation,  R2,  through  the  angle  i 
around  the  new  X  axis  (the  axis  through  the  ascending  node):  and,  (3)  a 
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rotation,  R3,  around  the  Z  axis  through  the  angle  Q.  This  gives 


fcostf 

-sin* 

O'! 

B12) 

Rl  3  Isintt 

cost? 

01 

l  0 

0 

1J 

r  1 

0 

0  1 

B13) 

po 

to 

II 

O 

cosi 

-sinil 

l  0 

sini 

cosiJ 

fcosQ 

-sinQ 

01 

B14) 

R3  =  IsinQ 

cosQ 

01 

l  0 

0 

1J 

If  R 

=  R3R2R1,  then 

Rll  3  cosQcosw  -  sinQcosisinw 

Rl2  =  -cosQsintf  -  sinflcosicosw 

Rl3  =  sinQsini 

R21  -  sinQcosw  ♦  cosQcosisinw 

R22  =  -sinQsinw  +  cosQcosicosw 

R23  *  -cosQsini 

R31  3  sinisinw 

R32  =  coswsini 

R33  =  cosi 


If  the  X-Y-Z  coordinate  system  is  the  true-of-date  geocentric 
inertial  system,  called  Celestial  Coordinates,  CC,  then  the 
transformation  to  Earth  fixed  coordinates  is  as  follows.  The  rotation 
of  the  Earth,  L,  since  the  reference  time  of  the  coordinate  system  is 

B15)  L  =  DFRADY*6 . 3003881  +  DELTDY*0 . 1720279  +  GHA 

where  DELTDY  is  the  number  of  days  since  the  reference  time  of 
the  coordinate  system, 

DFRADY  is  the  fraction  of  day 

GHA  is  the  Greenwich  hour  angle  at  the  reference  time 
of  the  coordinate  system. 
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II  the  Z-Y-Z  coordinate  rector  is  denoted  I  and  the  Earth-fixed 
coordinates,  X*,  Ye,  Ze,  are  denoted  X«,  then  X«  is  giren  by 


B16)  X*  =  TlI 


where 


f  cos(L  +  X)  sin(L  +  X)  O') 

B17)  Tl  =  I  -sin(L  +  X)  cos(L  +  X)  01 

10  0  1  J 


where  X  is  the  longitude. 

To  then  transfer  from  the  Earth  fixed  equatorial  system,  I«,  to  the 
topocentric  N-E-Y  (North,  East,  Vertical,  spherical  Earth,  geocentric) 
systea,  the  transformation  matrix  is  giren  by 


f  cos  (l  +  r/2)  0  sin(^  +  w/2)  1 

B18)  Tfc  =  I  0  1  0  1 

l-sin(l  +  r/2)  0  cos(^  +  r/2)  J 


where  9  is  the  geocentric  colatitude  and  ^  =  r/2  -  9  is  the  latitude. 
Combining  results,  we  then  hare 

f  cross') 

BIO)  X«  =  TlRI’  =  qi’  =  q’ I  down  I 

l  alongJ 

where  q  =  TlR  and,  in  the  notation  of  equations  B9)  and  BIO) ,  X*  is  the 
rector  (-down,  along,  -cross) ,  and  it  is  easy  to  show  that 

q’ll  *  -qi3  =  -cos(L  +  X)Ri3  -  sin(L  +  X)R23 

q*12  *  -qil  =  -cos(L  +  X)Rn  -  sin(L  +  X)R21 

q’l3  «  qi2  =  cos(L  +  X)Ri2  +  sin(L  +  X)R22 

0*21  *  -023  =  *in(L  +  X)Ri3  -  cos(L  +  X)R23 

q’22  =  -021  =  sin(L  +  X)Rn  -  cos(L  +  X)R21 

q’23  ®  022  ■  -sin(L  +  X)Ri2  +  cos(L  +  X)R22 

Q’31  =  -033  =  -R33 
Q’32  =  -031  =  -R31 
O’ 33  =  032  =  832 
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Comparison  of  BIO)  with  BIO)  shows  that  TGS  =  Q*  or  that 

Tx  =  Oil;  r7  *  021}  r*  =  031.  =  0l3J  ny  =  023>  n*  =  033 

and  that 

B20)  tanAn  =  ny/nx  =  Q23/Q13  *  cos(Q+L+A)/sin(Q-L-A) . 
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Effects  of  orbit  error  on  near-earth,  polar,  3 

satellite  aagnetic  data,  for  radial,  theta, 
and  phi  components  and  scalar  field  nagnitude. 

Temporal  distribution  of  obserratory  annual  9 

means  data  used  in  DMSP  models. 


Temporal  distribution  of  land  and  aeromagnetic  10 
data  used  in  DMSP  models. 


glO,  81^ »  and  hil  coefficients,  plus  Euler 
angle  and  bias  solutions,  generated  from  fits 
to  15  sets  of  DMSP  data  spanning  1984  to  1986. 
Table  4  data  is  plotted  in  Figure  7. 
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All  terms  for  MO,  MODM,  MODMS,  and 

MODMSDE (11 , 11 , 11)  models. 

Th^  relative  error  of  \n,  as  approximated  by 
equation  B6,  as  a  function  of  latitude  of  the 
DMSP  satellite. 
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Figure  Captions 


Figure  # 
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Caption 


Relation  of  Euler  angles  to  spacecraft  coordinate  system, 
describing  rotations  about  the  X,  Y,  and  Z  axes,  respectively. 
Rotations  around  X  and  Y  (figures  a  and  b)  are  left-handed, 
whereas  Z  is  right-handed. 

DMSP  orbital  X,  Y,  and  Z  magnetic  component  data  which  have 
had  an  estimated  field  model  removed,  revealing  strong 
periodicities  in  the  residuals.  The  dashed  line  is  a  spline 
fit  to  the  residuals. 

Power  spectra  of  X  and  Y  DMSP  residual  data  from  Figure  2. 

Y-coaponent  of  DMSP  data  from  Figure  2,  demonstrating  removal 
of  outliers,  magnetometer  rotation  and  bias  correction,  and 
subtraction  of  Fourier  periodic  function.  The  dashed  line  is  a 
spline  fit  to  the  residuals. 

World-wide  geographic  distribution  of  typical  three  day 
sequence  of  processed  DMSP  data.  Note  the  absence  of  data 
above  75°  latitude  and  the  large  mid-latitude  data  gaps  where 
the  torquing  coils  were  activated. 

Geomagnetic  field  spectra.  RQ  is  the  total  mean  square 
contribution  to  the  vector  field  from  all  harmonics  of  degree 
n.  The  spectra  are  from  field  models  produced  from 
uncorrected  and  corrected  DMSP  data;  the  Magsat  spectrum  is 
included  for  comparison. 


Plots  of  gi.0,  gil,  hjl,  Euler  angles,  and  biases  versus  time 
(yrs) ,  for  field  model  solutions  from  15  DMSP  data  sets 
spanning  1984  -  1986. 

Plots  of  ftiO,  fol,  and  fill  for  M0(ll,8,0),  110(11,11,0), 

MO  (11, 10,5)  and  MODM(11,10,5)  models,  at  epochs  1980  and 
1985.  Coefficients  derived  from  observatory  first 
differences  from  1940  through  1979  are  included  for 
comparison . 
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9  Plots  of  &0,  ftil,  and  fill  for  MO (11, 8, 5),  M0S(11,11,5) ,  and 
both  MOSDM  and  M0SDMDB(11,8,5)  and  (11,11,5)  models,  at 
epochs  1980  and  1985.  First  difference  coefficients  from 
1940  through  1985  are  included  for  comparison. 

10  Plots  of  glO,  gil  and  h].l  coefficients  for  individual  DMSP 
fits,  the  XGRF  1970,  1975  and  1980  models,  and  the  DE-2  1982 
calibration  model. 

11  Plots  of  the  gj.1  term  for  individual  DMSP  fits.  Figure  11A 
includes  the  DMSP  and  DE-2  calibration  models  for  comparison. 
Figure  11B  includes  the  M0DMSDE(11,9)  and  (11,9,5,3)  models, 
and  the  M0S(11,9)  model  for  comparison. 

B1  Relation  of  DMSP  orbit  with  local  n,  r,  v  coordinates  to 
earth-fixed  cartesian  coordinates. 

B2  Geometry  of  satellite  orbit  in  inertial  space.  X-Y-Z  is 

coordinate  system  with  Z  along  the  Earth’s  rotation  axis,  X 
along  the  vernal  equinox,  and  Y  so  as  to  give  a  right  handed 
system.  X’-Y’-Z’  is  the  spacecraft  coordinate  system  with  X’ 
radial,  V  normal  to  the  orbit,  and  Y’  so  as  to  give  a  right 
handed  system.  Y’  is  along  track.  The  inclination  is  i,  Q 
is  the  right  ascention  of  the  ascending  node,  and  u  is  the 
angular  position  of  the  spacecraft  measured  from  the 
ascending  node. 
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